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S^_ • ABSTRACT 

- ^ , We study Scliwinger-Dyson equation for fermions in Yukawa and Wess-Zumino models, 

X ■ 
^ • in terms of dynamical mass generation and the wavefunction renormalization function. 

In the Yukawa model with 75-type interaction between scalars and fermions, we find a 

critical coupling in the quenched approximation above which fermions acquire dynamical 

mass. This is shown to be true beyond the bare 3-point vertex approximation. In the 

Wess-Zumino model |]l|, there is a neat cancellation of terms leading to no dynamical 

mass for fermions. We comment on the conditions under which these results are general 

beyond the rainbow approximation and also on the ones under which supersymmetry is 

preserved and the scalars as well do not acquire mass. The results are in accordance with 

the non-renormalization theorem at least to order a in perturbation theory. In both the 

models, we also evaluate the wavefunction renormalization function, analytically in the 

neighbourhood of the critical coupling and numerically, away from it. 



1 Introduction 

Despite the success of quantum field theory in the description of the behaviour of ele- 
mentary particles in the perturbative regime of interactions, it still remains a challenge 
to understand the non-perturbative domain satisfactorily. One of the methods which has 
gained attention in this regard in recent years is the study of Schwinger-Dyson equations 
(SDEs) 0. Despite the difficulties involved in finding a non-perturbative truncation of 
these equations, this approach has been very successful in addressing issues like dynamical 
mass generation for fundamental fermions when they are involved in sufficiently strong 
interactions 0. Moreover, recent attempts, e.g., |^, |^, to improve the reliability of the 
approximations used have increased the credibility of the results obtained through such 
studies. 

Application of Schwinger-Dyson formalism to supersymmetric (SUSY) models has 
been less extensive. In supersymmetric Quantum Electrodynamics (SQED), based upon 
the arguments of non-renormalization theorem and gauge invariance, it is expected to be 
impossible to obtain dynamical mass generation for fermions [0 though some studies [§ 
argue that it is probably possible to break chiral symmetry dynamically in SQED. We 
postpone the study of SQED for a future work. In this paper, we take the simplest SUSY 
model, i.e the Wess-Zumino model and attempt to solve the corresponding Schwinger- 
Dyson equations for the fermion and scalar propagators. We believe this exercise will 
provide us with a deeper insight into how the role of supersymmetry in the context of 
dynamical mass generation translates into the language of Schwinger-Dyson equations. 

Such a study should provide us with a better starting point for more complicated 
SUSY theories such as SQED and SQCD. In the latter theory, a need also exists to further 
explore connections between the Holomorphic approach and that of the Schwinger-Dyson 
equations 0. 

We first study the Yukawa model with one real scalar and one Majorana fermion, 
which can be considered as a truncated Wess-Zumino model. We discuss this model in 
some detail for two reasons, first being that it is interesting in its own right because, 
after all, it is Yukawa interactions which are responsible for giving masses to fermions in 
the Standard Model (SM). Secondly, extending the Yukawa model by doubling the scalar 
degrees of freedom provides us with a clear understanding of how supersymmetry works. 



We use the quenched approximation. Keeping in mind the perturbative expansion of the 
3-point vertex beyond the lowest order and its transformation under charge conjugation 
symmetry, we propose an ansatz for the full effective vertex. One of the advantages of 
using this vertex is that the equations for the mass function M.{p'^) and the wavefunction 
renormalization F{p'^) decouple completely in the neighbourhood of the critical coupling, 
ac above which mass is generated for the fermions, and partly above it. We solve both the 
equations to find analytical expressions for -F(p^) and the anomalous mass dimensions in 
the neighbourhood of Uc- The results show that non-perturbative interaction of fermions 
with fundamental scalars can give masses to fermions in a dynamical way provided the 
interaction is strong enough. We use numerical calculation to draw Euclidean mass of the 
fermions as a function of the coupling, and confirm that it obeys Miransky scaling. We 
also evaluate F{p'^) numerically. We then extend the particle spectrum by doubling the 
number of scalars and imposing relations for the couplings that define the Wess-Zumino 
model. Due to the presence of the additional symmetry, we are able to extract useful 
information beyond the rainbow approximation. 

2 The Yukawa Model 

Consider a massless Lagrangian with one Majorana fermion and one real scalar interacting 
with each other through a 75-type interaction: 

C = ^{d,A)' + Iti^Yd,^) - \g'A' + tgi^^.M • (1) 

The corresponding Schwinger-Dyson equation for the fermion propagator, Sf{p), is dis- 
played in Fig. 1. Motivated by the success of the quenched approximation in QED 
and QCD, we neglect the fermion loops. Moreover, as a first step towards truncating 
the infinite set of Schwinger-Dyson equations, we drop all 4-point functions. The full 
scalar propagator can then be replaced by its bare counterpart. Using Feynman rules, the 
Schwinger-Dyson equation can be written as: 

-tS^\p) = -tS'f^p) - J-0^A-^91,) {^SAk)) {-2g^,TA{k,p)) (^) , (2) 

where q = k — p and Sp{p) can be expressed in terms of two Lorentz scalar functions, 
F{p'^), the wavefunction renormalization and A^(p^), the mass function, so that 



The bare propagator Sp{k) = l/(^), where the bare mass has been taken to be zero. We 
can project out equations for FijP') and A^(p^) by taking the trace of Eq. (2) having 
multiphed by ^ and 1 in turn. On Wick rotating to Euchdean space, 

where a = g'^/'iir. 

It is here that we cannot proceed any further without making an ansatzioi rA{k,p). 
Any ansatz for the 3-point vertex must fulfill at least the following requirements: 

• Perturbatively, we must have TA{k^p) = 1 + 0{g^). 

• It must be symmetric in k and p. 

Moreover, as the SDEs relate the 2-point function with the 3-point function, the expression 
for the full vertex is expected to involve functions FijP') or/and M.{p'^). 

The commonly used ansatz in the non-perturbative study of the SDEs is the bare 
vertex ansatz. For the Yukawa Model under discussion, it implies TA{k.,p) = 1. It agrees 
with the lowest order perturbation theory. The only truncation of the complete set of 
Schwinger-Dyson equations known so far that avoids any assumptions other than the 
smallness of the coupling at every level of this approximation is the perturbation theory. 
Therefore, it is natural to assume that physically meaningful solutions of the Schwinger- 
Dyson equations must agree with perturbative results in the weak coupling regime. It 
requires, e.g., that every non-perturbative ansatz chosen for the full vertex must reduce 
to its perturbative counterpart when the interactions are weak. Bare vertex fulfills this 
requirement to the lowest order in perturbation theory. Any other vertex which fulfills 
this condition and does not violate other requirements is at least as good as the bare 
vertex. One of the simplest non-perturbative vertices can be constructed by realizing 
that Eq.(^) yields the following expansion of F{p'^) in perturbation theory: 

Fi^p") = 1 + C(«) (6) 



Therefore, a simple candidate for the 3-point vertex can be written as: 



^A{k,p) 



F{k^)F{p'^) 



(7) 



Perturbatively, it gives 

^A{k,p) 



1 + C(«) 



[l + Oia)][l + Oia)] 

Therefore, to the lowest order in perturbation theory, our vertex ansatz reduces to the 
bare vertex. 

It is exceedingly complicated to ensure that at higher orders, the non-perturbative 
vertex reduces to its perturbative counterpart in the weak coupling regime. We do not 
aim at it in this paper. However, we demonstrate that even up to next to lowest order, 
i.e, to 0{a), our ansatz is correct to the extent that both the ansatz and the real vertex 
have the logarithmically divergent behaviour in the ultraviolet regime. 

Fig. 2 represents the perturbative expansion of the 3-point fermion-scalar vertex to 
0{a). Using Feynman rules, we can write it as follows: 



'^gin^A = -26-75 + 






(-25(75) iSpip - w) (-25(75) iSpik - w) (-2575) -^ (9) 



w^ 



which can be simplified to: 
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}t ^J(°) - (^7" + 7" fjj'^^ + K^""^ 
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(10) 

(11) 

(12) 
(13) 



{p — wY {k — wY 

The exact analytical expressions for these three integrals are known [|iy, |lT] . They involve 
basic functions of momenta k and p and a spence function. We believe that it is highly 
non-trivial to construct a non-perturbative vertex which reduces to this complicated form 
in the weak coupling regime. However, asymptotic behaviour can be reproduced to some 
extent. Simple power counting reveals that the integrals J^^-* and JJ^' are perfectly well- 
behaved in the ultraviolet regime. However, K'^'^'> is logarithmically divergent. We now 



show that our vertex ansatz also exhibits this behaviour. Using the fact that perturba- 
tively M.{p'^) = 0, we can re- write Eq. (4) as follows: 






(14) 



where we have employed Eq. (6), and the Feynman rule for the vertex. Carrying out 
angular and radial integrations respectively and retaining the leading log terms, we get 

n2 



F{p' 



1 + — n 4^ 

27r A2 



(15) 



Therefore, the proposed vertex ansatz can be written perturbatively as follows 

1 



^A{k,p) 



l + -ln^ + 0{a^) 



(16) 



F{k'^)F{p^) ^ ' tt'" A2 

which is logarithmically divergent in the ultraviolet regime just as the real vertex to 0{a). 
Therefore, perturbatively our vertex ansatz is more realistic than the bare vertex. 

An added advantage of using the proposed vertex ansatz is that Eq. (^ can be solved 
independently of Eq.(^. Therefore, it serves a purpose similar to that of Mandlestam's 
choice [|1^] of the 3-gluon vertex in studying the Schwinger- Dyson equation of the gluon 
propagator. 

As the unknown functions F and A4 do not depend upon the angle between k and p, 
we can perform angular integration to arrive at 



F{p') 
M{p') 



1 
a 

TT 



1 



^ fdk^ 

27tJ k^ + M^{k^) 



k' 



-Mp' - k^) + 0{k^ - p") 

p4 



dk^ 



M{k^ 



pz 



9{p^ - e) + 9{k' 



p 



(17) 
(18) 



k^ + M^{k^ 

Such equations are known to have a non-trivial solution for the mass function above a 
critical value of the coupling a = ac- In the neighbourhood of the critical coupling, when 
the generated mass is still small, we can put Ai^ = 0. Then Eqs. (|T^Jl8| ) decouple from 
each other completely. The leading log solution for F{p^) is then: 



F(/ 



27r A2 



(19) 



As for the mass function, multiplicative renormalizability demands a solution of the type 
M(p2) ~ (p2)-'5. Substituting this in Eq. (|18|) and performing radial integration, we find 
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1 1 / a , , 

where ac = tt/A. For a > ac the solution of the mass function enters the complex plane 
indicating that a phase transition has taken place from perturbative to non-perturbative 
solution corresponding to the dynamical generation of mass. Numerically, above ac, we 
solve Eq. ([T8|) in a two-step process. We first use the iterative method to get close to 
the solution and then refine the answer by converting the integral equation into a set of 
simultaneous nonlinear equations to be solved by Newton- Raphson method. In Fig. 3., 
we have drawn the Euclidean mass M(which can be taken to be A^(0)) as a function of 
the coupling a. We see that it obeys Miransky scaling law and can be fitted to the form 



M 

- = exp 



^^B 



(21) 



very well by the choice A = 0.977r and B = 1.45. These numbers are close to the ones 
found in [l^ although the value of the critical coupling is of course different. The slight 



mismatch in the values of A and B is due to the fact that in the logarithmic grid of 
momenta, we choose 30 points per decade and do not extrapolate the result to an infinite 
number of points, an exercise carried out in [|13|. We also compute -F(p^) for various 
values of a and find that the closer we approach ac, where the generated mass is still 
small, starting from a larger value of a, the numerical result gets closer and closer to the 
analytical result as expected. Fig. 4. 

3 The Wess-Zumino Model 

We now extend the particle spectrum by doubling the number of scalars to discuss the 
massless Wess-Zumino model, characterized by the following Lagrangian: 

£ = ^{d.Af + ^{d.By + ^{^^Pr^,i;) - \g\A' + B'f - gi,{B - lA^,)^, . (22) 

The Schwinger- Dyson equation for the fermion propagator in this model is depicted in 
Fig. 5. Before we embark on solving this equation, we must re-address the validity of the 
ansatz for the 3-point vertex which we made for the Yukawa case. To the one loop level 
in perturbation theory, one now has contributions from both the scalars as depicted in 



Fig. 6. Therefore, e.g., for scalar A, we can write 

/d tv i 

7^^ (-2^75) iSpip - w) (-2^75) iSpik - w) (-2^75) — 

/nil] j 

The same is the case for Tb, i.e., in the presence of both the scalars A and B, none 
of the 3-point vertices gets modified at 0{a) in perturbation theory. Therefore, in the 
Wess-Zumino case, it is more reasonable to use the bare vertex instead of the ansatz 
earher made. Though one would now expect to solve coupled integral equations for Fijp') 
and A4(p^), a miraculous cancellation of terms takes place as evident from the following 
Schwinger-Dyson equation for the fermion propagator: 



iFip^) i TT^ 



F(fc^) 1 



+ 4 I d^k 



F{e) i_ 



(24) 



Taking the trace of this equation, we get 

M{p^) = 

As the cancellation of terms takes place at the very beginning, it is easy to see that 
dynamical mass generation will remain an impossibility for the full vertex and the full 
scalar propagator as long as they are identical for both the scalars. The vertex corrections 
for A and B have been proven to be equal up to O(a^) |T^. We shall shortly see that the 
same is true for the full scalar propagator at least up to 0{a). This is in accordance with 
the arguments based on non-renormalization theorem. SUSY plays a role in providing 
same number of bosonic and fermionic degrees of freedom. We have seen from the case of 
the Yukawa Model that without this equality, it will not be possible to prevent dynamical 
mass generation. Secondly, SUSY imposes relations on couplings of the two scalars with 
the fermions. This relationship is crucial in preventing dynamical mass generation. 
As far as wavefunction renormalization F{p'^) is concerned, its leading log behaviour gets 
modified slightly, by the inclusion of the other scalar, to: 

F(/) = l + -ln^ . (25) 



Although it is an interesting conclusion in its own right that supersymmetry prevents 
dynamical mass generation for fermions in the Wess-Zumino model, another important 
issue to probe will be whether supersymmetry itself remains intact, i.e., whether the 
scalars can also be kept massless. This is what we discuss now. The Schwinger-Dyson 
equation for the scalar (for example A) has been depicted in Fig. 7. A scalar propagator, 
unlike a fermion, needs only one unknown function to describe it. But we shall prefer to 
split it into two parts and write the full scalar propagator as follows: 



SAip) = .. TJ,., ■ (26) 



FAf 
p^ — M^ip"^ 

The non-zero value of the mass function A^^(p^) will be responsible for shifting the pole 
from p^ = to some finite value, generating the mass for the scalar dynamically. FaIp"^) 
on the other hand is the scalar wavefunction renormalization. The SD-equation for the 
scalar propagator in Euclidean space can now be written as: 



F4(p2) ^ ' 2n3j k^ + Ml{k^) 271^ J k^ + Ml{k^) 

2a r,^ J{k^)F {q^ 
vr^ J k'^q 



^ d'k ',i^;'' ^A{k,p)k.q (27) 



where we have used the fact that the fermions do not acquire mass. If we want to preserve 
supersymmetry and do not want the scalars to acquire mass, we must have: 



MAip')=MBip') = . (2^ 



We are then left with: 



a f d k 



Fa{p^) 27r^p^J A;2 

and there is a similar equation for the scalar B: 



SFAik') + Fsik') - 4^ F{e)F{q')TA{k,p) 

q^ 



(29) 



1 a r d^k 

Fb{p^) ^ ^ 27r3p2 J ^ 



SFnik') + FAik') - 4^ Fik')Fiq')Tnik,p) 

q- 



.(30) 



These equations should yield a solution for Fa{p'^) and Fb{p'^) such that it does not change 
the position of the pole for the scalar propagator and that the quadratic divergences cancel. 
It is well-known that it does happen in perturbation theory to 0{a). In fact, one can 
evaluate Fa{p'^) and Fb{p'^). The leading log expression for these functions to 0{a) is 

9 



FA{p')=Fsip') = l + -\n^ (31) 

which is exactly the same expression as that for F{p'^) for the fermion propagator. This 
result indicates that supersymmetry need not be broken. 

4 Conclusions 

We have studied the Schwinger- Dyson equations for the Yukawa (a scalar interacting 
with a fermion with a 75 type interaction) and the Wess-Zumino models. In the simple 
Yukawa model, we propose a vertex ansatz which we argue should perform better than 
the bare vertex. In the quenched approximation, we find dynamical mass generation for 
fermions above a critical value of the coupling Uc = 7r/4. The generated Euclidean mass 
obeys Miransky scaling. When we extend this Yukawa model to equate the scalar and 
fermionic degrees of freedom (Wess-Zumino model), we find that a neat cancellation of 
terms occurs and there is no mass generation for the fermions. This fact remains true 
beyond the rainbow approximation and is supported by perturbative calculations available 
for the 3-point vertex to 0{a'^) and of the scalar propagator. This result was expected on 
the basis of non-renormalization theorem. The two approaches will remain in agreement 
provided the full vertex and the full scalar propagator as long as they are identical for 
both the scalars to higher orders in perturbation theory as well. 

If supersymmetry has to be preserved, the scalars should also acquire no mass dynam- 
ically. Studying the Schwinger-Dyson equations for the scalars, we observe that such a 
solution is allowed and in fact leads to the wavefunction renormalization function for the 
scalars which is exactly the same as that for the fermion. 

It is more interesting to see the role of supersymmetry in more complicated theories 
such as SQED. The studies so far carried out in superfield and component formalism seem 
to arrive at different conclusions. We plan to present our work in this context in a future 
publication. 
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Figure Captions 

Fig. 1. Schwinger-Dyson equation for the fermion propagator in the Yukawa model. The 
solid lines represent fermions and the dashed the scalar. The solid dots indicate full, as 
opposed to bare, quantities. 

Fig. 2. One loop perturbative expansion for the vertex in the Yukawa model. 

Fig. 3. The dynamically generated mass (M/A) versus the 3-point coupling a in the 
Yukawa model. The critical coupling is a^ = 7r/4, above which the mass can be seen to 
be bifurcating away from the chirally symmetric solution, os represent the numerical result 
and +s the numerical fit to the form M = Aexp —A/Ja/ac — 1 + B with A = 0.97n and 
B = 1.45. 

Fig. 4. The wavefunction renormalization function F{p'^) in the Yukawa model for various 
values of the coupling a. The solid line corresponds to the analytical expression 1 — a/An + 
{a/27r) ln(p^/A^) in case of no mass generation for a = 0.78. 

Fig. 5. Schwinger-Dyson equation for the fermion propagator in the Wess-Zumino model. 

Fig. 6. One loop perturbative expansion for the vertex in the Wess-Zumino model. We 
have shown the case for scalar A. A similar diagram exists for scalar B. 

Fig. 7. Schwinger-Dyson equation for the scalar propagator in the Wess-Zumino model. 
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